Three Orbital Model for the iron-based superconductors 
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The theoretical need to study the properties of the Fe-based high-Tc superconductors using reliable 
many-body techniques has highlighted the importance of determining what is the minimum number 
of orbital degrees of freedom that will capture the physics of these materials. While the shape of 
the Fermi surface (FS) obtained with the local density approximation (LDA) can be reproduced by 
a two-orbital model, it has been argued that the bands that cross the chemical potential result from 
the strong hybridization of three of the Fe 3d orbitals. For this reason, a three-orbital Hamiltonian 
for LaOFeAs obtained with the Slater-Koster formalism by considering the hybridization of the As p 
orbitals with the Fe d^z, dyz, and d^y orbitals is discussed here. This model reproduces qualitatively 
the FS shape and orbital composition obtained by LDA calculations for undoped LaOFeAs when 
four electrons per Fe are considered. Within a mean-field approximation, its magnetic and orbital 
properties in the undoped case are here described for intermediate values of J/U . Increasing the 
Coulomb repulsion U at zero temperature, four different regimes are obtained: (1) paramagnetic, 
(2) magnetic (vr, 0) spin order, (3) the same (tt, 0) spin order but now including orbital order, and 
finally & (4) magnetic and orbital ordered insulator. The spin-singlet pairing operators allowed 
by the lattice and orbital symmetries are also constructed. It is found that for pairs of electrons 
involving up to diagonal nearest-neighbors sites, the only fully gapped and purely intraband spin- 
singlet pairing operator is given by A(k) = /(k) ^^^^ '^k,a,t'^-k,a,4, with /(k) = 1 or cos fcj, cos fc^ 
which would arise only if the electrons in all different orbitals couple with equal strength to the 
source of pairing. 

PACS numbers: 71.10.-w, 71.10.Fd, 74.20.Rp 



I. INTRODUCTION 

The discovery of high-Tc superconductivity in a fam- 
ily of iron-based compoundsir— is offering a new concep- 
tual framework to study the non-standard pairing mech- 
anism that seems to induce these exotic superconducting 
statesi^ The magnetism present in several of the parent 
compoundsiSrJd and the high critical temperatures are 
reminiscent of properties observed in the cupratesJ^ But 
there are also clear differences, since the parent com- 
pounds are metallioiSrJi and the Fermi surface (FS) is 
determined by more than one orbitalJ^"— The multior- 
bital nature of the problem poses a challenge to the de- 
sign of minimal models that can be studied with powerful 
techniques, such as numerical methods. Ab-initio calcu- 
lations making use of the local density approximation 
(LDA) indicate that the five id orbitals of Fe strongly 
hybridize to form the bands that are close to the chem- 
ical potentialfi^"— However, the FS appears to be de- 
termined by bands that have mostly dxz and dyz char- 
acter, and this observation has been confirmed by polar- 
ized angle resolved photoemission spectroscopy (ARPES) 
experimentsi^ This supports the notion that the mini- 
mum number of orbitals to be considered to study the 
pnictides could be two. In fact, a two-orbital mini- 
mal model has been proposed22i2^ and it has been stud- 
ied with numerical techniques on small clusters ^^i^ as 
well as with several other approximations , ^5- _30 ggth 
numerica P^i^^ and mean-field^^ calculations indicate that 
the magnetic metallic regime observed experimentally in 
the undoped compoundsiSril is stabilized for intermedi- 



ate values of the Coulomb repulsion U and the numerical 
calculations suggest that, upon doping in this regime, 
the most favored pairing operator is interorbital and has 
symmetry B2gM^ 

On the other hand, several authors have claimed that 
a two-orbital model may miss important features of the 
real systemJ^ii^"— It has been argued that a mini- 
mal model for the pnictides should contain at least three 
orbitals for mainly two reasons: (i) relatively small 
portion of the electron-pocket FS of LaOFeAs is deter- 
mined by a band of mostly d^y character and (ii) the 
bands that produce the two hole pockets should be de- 
generate at the center of the Brillouin zone (BZ), which 
is not the case when only two orbitals are considered. 
The important question is how much these shortcomings 
of the model impact the most relevant properties of the 
pnictides. The aim of this paper is to construct a three- 
orbital model that addresses these concerns and compare 
its properties with the two-orbital case. This is important 
because in other areas of condensed matter physics, such 
as the manganites, we have learned that a simple single- 
orbital model is often sufficient to capture qualitatively 
the phenomenon of colossal magnetoresistancC)^ while 
clearly a two-orbital model is still necessary to properly 
describe additional properties such as the magnetic and 
orbital order observed in these materialsi^ Similarly the 
minimal two-orbital model for the pnictides appears to 
reproduce the experimentally observed magnetic order, 
and it is interesting to investigate to what extent the 
inclusion of the third orbital xy modifies these results. 

While the pnictides exhibit a structural phase tran- 
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sition at similar or slightly higher temperature as the 
onset of antiferromagnetic (AF) order, experiments have 
not yet addressed the issue of orbital order and also have 
not provided consensus regarding the symmetry of the 
pairing operatori^"— An investigation of the magnetic 
and orbital orders, as well as pairing symmetries, that 
arc allowed in a three-orbital model compared to those 
that arc possible in the two-orbital case will shed light on 
the importance of the role that the additional d^y orbital 
should play in theoretical discussions. 

This paper is organized as follows. In SecHIl the three- 
orbital model is introduced. Section Hill contains a mean- 
field analysis of the magnetic and orbital properties of 
the undoped system. Section |IV] is devoted to the classi- 
fication and analysis of the spin-singlet pairing operators 
allowed by the orbital and lattice symmetries. A sum- 
mary and conclusions are presented in Sec. [V] 



II. THREE ORBITAL MODEL FOR THE 
PNICTIDES 

As explained in the Introduction, it has been suggested 
that at least three orbitals may be needed to describe the 
superconducting pnictides because the bands that deter- 
mine the Fermi surface of LaOFeAs are mostly composed 
by orbitals d^z, dyz, and dxy'^^^^—'^^^ The notation 
xz, yz, and xy will be used for these orbitals, respec- 
tively, for better readability. The need to include at least 
three orbitals in a realistic model was first pointed out 
in Ref. |3l| where a three-orbital model was constructed 
using the symmetry properties of the Fe-As planes and 
LDA results. A shortcoming of that proposed model was 
that it contained an spurious hole-pocket FS around the 
M point in the extended BZ notation. It was argued^ 
that a fourth orbital should be added to remove the extra 
pocket. However, it will be shown in Sec. Ill Bl that this 
spurious pocket can actually be removed entirely within 
a three-orbital model formalism, i.e. without adding a 
fourth orbital. 

One important issue that needs to be addressed is the 
electronic filling to be used in a three-orbital model. 
Band calculations have determined that the undoped 
pnictides contain six electrons distributed among the five 
3d orbitals of each Fe atom. One procedure to determine 
the filling for a model with a reduced number of orbitals 
is to start from the crystal field splitting and fill the levels 
accordingly from the lowest energy up. In the two-orbital 
model that considers only the xz and yz orbitals, such a 
consideration would suggest that half-filling is the correct 
electronic density;^"— because the — y^ and 3z^ — 
orbitals are assumed to be fully occupied with four of the 
six electrons and the xy orbital is assumed empty, leav- 
ing two electrons to populate the xz and yz orbitals. In 
addition, this filling assignment is the only one that al- 
lows to reproduce the LDA calculated FS. Applying the 
crystal-field splitting rationale to the three-orbital model 
with xz, yz, and xy orbitals, this argument leads to a 



filling of one third (i.e. two electrons in three orbitals) <^ 
However, for such a filling we have not been able to re- 
produce the LDA shape of the FS. Thus, the filling must 
be adjusted to approximately reproduce the FS and the 
orbital occupation numbers obtained with LDA. In fact, 
band structure calculations suggest that the three-orbital 
system should be more than half-filled and actually have 
a filling of roughly two thirds (i.e. four electrons in the 
three orbitals) j^i^ Our analysis shows that a FS with ap- 
proximately a similar size for the hole and electron pock- 
ets can be obtained both at fillings around one and two 
thirds (i.e. two and four electrons in the three orbitals), 
but the two almost degenerate hole-pockets around F de- 
mand a filling larger than half-filling. Thus, the focus of 
our effort will be on a filling of 2/3, as in Ref. [s^. As 
it will be discussed below, non-half-filled orbitals allow 
the orbital degree of freedom to be active and actually 
it has been argued that orbital ordering phenomena may 
play a role in these materialsi^ This orbital order is un- 
likely in the half-filled case which is the natural filling for 
the two22i2^ and four— orbital models for the pnictides. 
Once again, note that some authors have considered half- 
filling in the three-orbital case?^ but this leads to the 
"unwanted" hole pocket around M. 



A. Real Space 



To construct the tight-binding portion of the three- 
orbital Hamiltonian for the pnictides, the Slater-Koster 
procedure described in Ref. will be followed. Nearest- 
neighbor and diagonal next-nearest-neighbor hoppings 
will be considered for all the orbitals. It is clear that the 
hopping terms for the xz and yz orbitals arc the same as 
in the previously discussed two-orbital model. 



^rrz.yz ^ _^ y-/^t ^ ^^t ^ +h.C.) 

i,cr 

- 12 y^(d^- d-^- + £ d.^^ + h.c.) 

i,fT 

(d^- d.^^^^ -\-d^. -\-h.c.) 

/ J \ \,xz,a i-\-fi-\-v,xz,a ' i,yz,(T i-\-fi-\-i',yz,a ' / 

1,(1,0,(7 

''^^^^l.xz.a^i-\-x-\-y,yz,a ^" ^l,yz,cr^i-\-x-\-y,xz,(7 ^■^•) 
i.cr 

- tAy^id^. d.^^ ^ + d.^^ ^ + h.c.) 

^ / i.xz.a i-i-x — y.yz.a ' i,yz,(T i-^-x — y.xz.a ' / 
i,fT 

- fJ.^{ni^xz + ni,yz), 

(1) 
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while the intra-orbital hoppings for the xy orbital are 
given by 



J2 i4.xy.A+il+0.xy.a + '^-' 



(2) 



where lS.xy is the energy difference between the xy and 
the degenerate xz/yz orbitals. The hybridization be- 
tween the xz/yz and the xy orbitals is given by 



jjxz,yz;xy 



i.a 

^^E 

i.a 

^^E 
^sE 

i,(T 

i,(T 
i,(T 
i.a 

^sE 

i,(T 

isE 

i,(T 

tsE 



(-1 
(-1 
(-1 
(-1 
(-1 
(-1 
(-1 
(-1 
(-1 
(-1 
(-1 
(-1 



i^xz^a i--\-x,xy.(T 



d}. d. ^ 

i,xy,a i'\'X,xz,(j 



^i^yz,(7^i~\-y,xy,a 



+ h.c] 
+ h.c] 
+ h.c] 
+ h.c] 



d 

i,xz,(T i-\-x-\-y,xy,a 



7t 



'i,xy,cr i~\-x~\-y,xz,(7 



i, a; 2,(7 y,a;y,cr 



:I d.^~ - 



i,yjs,f7 i-)-i+y,a:;y,(T 



i.xy .a i-^-x-^-y .y z 
1 , y ^ , (T 1 + a; — y , a: y , cr 



/i.e.] 
/i.e.] 
/i.e.] 

f /i.e.; 

/i.e.] 
/i.e.] 
/i.e.] 



i.xy .(7 \-\'X~y,yz 



„ + h.c]. 



(3) 



The hopping parameters ti in Eqs.(IT]l3]) will be deter- 
mined by fitting the band dispersion to band structure 
calculation results. The chemical potential ^ is set to 
a two-thirds filling, as already discussed. The opera- 
tor d\^ ^ (d; ^ ^) creates (annihilates) an electron at site 
i, orbital a — xz,yz, xy, and with spin projection a. 
Tii.a = n-i^a.-]- + n-i^a.i SLve the corresponding density op- 
erators. Previously proposed three-orbital models^i^ 
only contained the nearest-neighbor (NN) hybridization 
tj, but since next-nearest neighbor (NNN) terms are in- 
cluded for the intra-orbital component, as well as for the 
hybridization between xz and yz, they should also be in- 
cluded in the hybridization with xy. In Sec. Ill Bl it is 



TABLE I: Parameters for the tight-binding portion of the 
three-orbital model Eqs.((S]) to pO|) . The overall energy unit 
is electron volts. 



ti 



t2 h 



0.02 0.06 0.03 -0.01 0.2 0.3 -0.2 -17/2 0.4 



shown that these NNN terms with hopping /§ turn out 
to be crucial to provide the proper orbital character for 
the electron pockets when compared with LDA results. 
Finally note that the hybridization Eq. ([3]) contains fac- 
tors (—1)''' that arise from the two-iron unit cell of the 
original FeAs planes. 



B. Momentum Space 



The Hamiltonian Eqs. can be transformed to mo- 



mentum space using ,^ = ^ e '^''^d\,a,a^ where 
k is the wavevector and N the number of sites. Note that 
the Fourier transformed Hamiltonian is defined in the ex- 
tended or unfolded BZ.-^i^ As pointed out in Ref. [sil, 
the real-space Hamiltonian presented in Eqs. is in- 

variant under a translation along the x or y directions 
followed by a reflexion about the x-z plane. When the 
eigenstates of the Hamiltonian are labeled in terms of 
the eigenvalues of these symmetry operations, then the 
momentum-space Hamiltonian can be expressed in terms 
of a pseudocrystal momentum that will be called k that 
expands the unfolded BZ that corresponds to a single 
Fc unit-cell in real space. Note that folding the three 
bands used here into the reduced diamond-like cell de- 
fined in the first quadrant hy kx + ky < tt, and symmetric 
points in the other three quadrants, doubles the number 
of bands to six as expected. 

Equations (IT][3]) become in momentum space 



i7TB(k) = E T^''ik)dl^,J^ 



(4) 



with 



T = 2t2 cos kx + 2ti cos ky + At^ cos kx cos ky — fj.(5) 
T^^ ~ 2ti cos kx + 2t2 cos ky 4- 4t3 cos kx cos ky — ^{6) 

^ 2t5(cOS/Ca: -f COsfcy) 

-1-4^6 COS /Ca; cos /Cj, — /i -|- Aa:j,, (7) 

J.12 ^ T^i 4i4sinfc^sin/i:y, (8) 
= T^^ = 2it7 sin + Aits sin kx cos ky, (9) 
^ T^^2itr sin ky + Aits sin ky cos kx , (10) 

where a bar on top of a matrix element denotes the com- 
plex conjugate. Since the Hamiltonian for a one-iron unit 
cell has been considered, then k runs within the corre- 
sponding extended BZ — tt < kx,ky < tt. 

The hopping parameters /, and Axy were chosen to 
reproduce the shape of the LaOFcAs FS obtained using 
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FIG. 1: (Color online) (a) Band structure and (b) Fermi sur- 
face of the tight-binding (i.e. non- interacting) three-orbital 
model, with parameters from Tab. |T] and in the unfolded BZ. 
The diagonal thin solid line in |(b)| indicates the boundary of 
the folded BZ. In panels (c-e), the orbital contributions to 
the two hole and one of the electron pockets are given. The 
winding angle 9 is measured with respect to the fcy-axis. The 
second electron pocket is analogous to the one shown simply 
replacing xz by yz. In all panels, the dashed lines refer to the 
xz orbital, the solid to yz, and the dotted to xy. 



LDA calculations. These parameters are given in Tab. ID 
The chemical potential ^ = 0.212 was fixed to this value 
to ensure a filling of two thirds, as discussed in Se c. Illl 
The resulting band dispersion is presented in Fig 



1(a) 



along high symmetry directions in the extended BZ and 
the corresponding FS is shown in Fig. 1(b) The two 



hole-pocket FSs, labeled ai and a2, are formed by two 
bands that are degenerate at F in agreement with LDA, 
and both of them are found around the F-point instead of 
M in the extended zone. Thus, one of the shortcomings 
of the two-orbital model has been corrected. It can also 
be observed that there is no hole-pocket FS around M 
which was a problem encountered in Refs. [siUs^ 

The orbital composition of the hole pockets is displayed 
in Figs. [1] (c) and (d). They are given by a linear com- 
bination of the xz and yz orbitals, in agreement with 
LDA^s^ and ARPES resuhs,^ The electron pockets at 
X {Y) arise from a linear combination of the yz [xz) 
and xy orbitals. The orbital composition for the electron 
pocket around X, labeled /3i in Fig. |l(b)[ is displayed 



in Fig. 1(e) here it can be observed that the orbital 



character changes from purely xy along T-X to predom- 



inantly yz along X-M , as predicted by LDA^i^ and also 
found with ARPES techniques!^ As can be deduced from 
Eqs. ^ and pU]) . setting tg, = —t-ijl ensures that the 
electron pockets have pure xy character along the V-X 
and Y-Y directionsi^ A large yz (xz) contribution along 
X-M (F-A/) requires the hybridizations t-j and to be 
quite robust)^ 

In summary, a tight-binding Hamiltonian has been 
constructed that captures the generic shape and orbital 
composition of the FS for undoped LaOFeAs by consider- 
ing only the three xz, yz, and xy orbitals, and assuming 
that they share four of the six electrons that populate 
the five Fe Sd levels. 

The Coulombic interacting portion of the Hamiltonian 
is given by: 

i,Q i,a</3 

— 2 J Si_Q, • Si^^ 

i,a<)9 

i,a</3 

(11) 

where a,/9 = xz,yz,xy denote the orbital, Si_Q is 
the spin (electronic density) in orbital a at site i, and 
the relation U' = [/ — 2 J between these Kanamori pa- 
rameters has been used (for a discussion in the manganite 
context see Ref . [32] and references therein) . The first line 
terms give the energy cost of two electrons located in 
the same orbital or in different orbitals on the same site, 
respectively. The second line contains the Hund's rule 
coupling that favors the ferromagnetic (FM) alignment 
of the spins in different orbitals at the same lattice site. 
The "pair-hopping" term is in the third line and its cou- 
pling is equal to J by symmetry. The values used for U 
and J can be substantially smaller than the atomic ones, 
because the interactions may be screened by bands not 
included in our Hamiltonian. These Coulombic interac- 
tion terms have been used and discussed in several pre- 
vious publications2^"— where more details can be found 
by the readers. All the energies in this paper are given 
in electron volts. 



III. NUMERICAL RESULTS 

A. Mean-field Approximation and Ordered Phases 

To investigate the magnetic and orbital order prop- 
erties of the three-orbital model in the presence of the 
onsite Coulomb repulsion, here a mean-field approxima- 
tion will be used. As it was carried out for the two- and 
four-orbital models in Ref. [2^, here only the mean-field 
values for diagonal operators will be considered^ (for de- 
tails see App. E| . In the present case, the degenerate xz 
and yz orbitals make up most of the weight close to the 
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TABLE II: Magnetic and orbital ordering wavevectors for the 
possible ordered phases discussed in the text. The third col- 
umn indicates the panel of Figs. [2]and[3l where a schematic 
sketch for the corresponding ordered phase can be found. 

panel in Figs. [2] or |3] 
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FIG. 2: (Color online) Cartoons representing the spin- and 
orbital-order configurations considered in this effort. Since 
four electrons occupy three orbitals, perfectly ordered states 
have one doubly occupied orbital and the unpaired electrons 
are located in the other two orbitals, each singly occupied, 
forming a net spin S = 1. In this cartoon, the orbital drawn 
indicates the doubly occupied one and the arrows indicate the 
orientation of the total spin at that site. |(a)| (tt, 0)-spin and 
ferro-orbital (FO) order favoring the y z orb ital, |(b)| (7r, 0)-spin 
and FO order favoring (xz + yz)/\/2, |(c)| (tt, 7r)-spin and FO 
order. (d)|(f) states with the sam e ordering vector for the 
spins and the orbitals: (7r,0) for |(d)| and (e) (7r,7r) in (f) (g) 
(tt, 0) for spins an d (0 , tt) for orbitals; (hj "(tt, 0) for spins and 
(tt, n) f or or bitals; (i) {n,n) for spins and {tv,0) for orbitals; 
(b) and | (e) | illustrate phases where the orbital order does not 
feature alternating xz and yz orbitals, but the combinations 
{xz + yz)/y/2 and {xz — yz)/y/2 do. 



(a) 



(b) 



FIG. 3: Color online) Cartoons for the spin ferromagnetic 
phases taken into consideration in this effort. Results are 
depicted with the convention of Fig. [21 but for ferromagnetic 
spins combined with (a) FO, (b) (7r,0), and (c) (tt, tt) orbital 
order. 



chemical potential. In previous investigations of the two- 
orbital xz-yz model, very weak orbital order has been 
found within a mean-field approximation at intermedi- 
ate Coulomb repulsion in the half-filled case, and none in 
the large-?/ limit.~ Numerical simulations did not pro- 
vide indications of orbital order at half-filling eitheri^ 
However, the xz and yz orbitals are not expected to be 
half filled in the present case: In the noninteracting t/ = 
limit, the xy orbital is approximately half filled, so that 
the subsystem consisting of the two degenerate xz and yz 
orbitals turns out to be approximately 3/4 filled. Con- 
sequently, orbital order within the xz / yz-subspace could 
now occur, analogously to the case of quarter filling in 
the xz/yz subsystemi^ 

Three possible orbital-order patterns will be consid- 
ered: (i) Ferro-orbital (FO) order which corresponds to 
the orbitals xz and yz having different electronic densi- 
ties, (ii) alternating orbital (AO) order, and (in) stripe 
orbital (SO) order. Combined with the magnetic spin 
order, these orbital orders lead to a large variety of pos- 
sible combinations of polarized or alternating spin and 
orbital order Here, phases that can be expressed using 
(at most) two ordering vectors have been considered, i.e. 
qi for magnetic order and q2 for orbital order. The ex- 
pectation values for the mean-field proposed states can 
then be expressed as 



71 



r,xy,a/ 



a 

( 

2 



•'xy 



(12) 
(13) 



where the first equation with the mean-field parameters 
nxyi^xy describes the xy orbital and the second equa- 
tion with parameters n, m,p, and q applies to the xz/yz 
subsystem, with a = ±1 indicating the xz/yz orbitals. 

The fact that the xz/yz space is not SU{2) sym- 
metric introduces another degree of freedom in addition 
to the ordering vector: Ferro-orbital order (i.e., site- 
independent orbital densities) can favor either the xz or 
yz orbitals [as shown in Fig. 2(a)| , or any linear combina- 
tion 10) = cos (f>\xz) + sin (j)\yz). As an example, Fig. 2(b) 
illustrates a state with </ 
metric combination (\xz 



7r/4 corresponding to a sym- 
\yz))/ \/2. The same holds for 
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alternating orbital order: alternating order correspond- 
ing to 4> = 7r/4 can be seen in Fig. 2(e) Consequently, 



mean-field calculations were performed for several values 
of for each phase. The xy orbital might be similarly in- 
volved in such a linear combination, because the crystal- 
field splitting separating it from the xz and yz orbitals is 
not very large. However, Exact Diagonalization of 2 x 2 
clusters did not give any indication for such behavior: 
the dominant states in the low-energy eigenstates that 
we analyzed involved almost exclusively singly occupied 
or unoccupied xy orbitals. Moreover, the xy orbital is al- 
most singly occupied at C/ = and has low weight at the 
Fermi energy, so that ordering phenomena in the most 
relevant intermediate- C/ regime may be expected to in- 
volve mostly xz and yz orbitals. The states considered 
are shown in Figs. [5] and [3] and the corresponding values 
of (\i are given in Tab. |lll Which spin-orbital pattern is 
stabilized depends on the interaction parameters U and 
J, as described in Sec. IIII Bl The phase with the largest 
stability range, and which is stable for the most realis- 
tic parameter choices, is the (tt, 0)-AF orbital disordered 
(OD) phase, which will be discussed in more detail in 
Sec HlTC] 



B. Magnetic and Orbital Orders in the Undoped 
Regime 

Magnetic order with wavevector qi = (tt, 0) (or (0, tt)) 
and OD was found to be stable in a broad - and especially 
the most realistic - range of interaction parameters. How- 
ever, the phase diagram in the J /U vs. U plane turns out 
to contain a large variety of metallic disordered phases, 
metallic phases with different kinds of magnetic and/or 
orbital order, and insulating magnetically and orbitally 
ordered regions. Figure S] shows a qualitative rendition of 
the resulting phase diagram. A more accurate quantita- 
tive determination of the boundaries as well as a detailed 
description of all phases will be presented in a future pub- 
lication. For a realistic Hund's rule coupling J = f7/4, we 
discuss the properties of the ground state in the different 
regimes encountered varying U in Sec. IHI Cl below. 

It has been found that at large interaction strengths U , 
after the magnetic order with qi = (tt, 0) is established, 
alternating- [q2 = (tt, tt)] and ferro-orbital [q2 = (0, 0)] 
order develops for 1/5 < J/U < 1/3. For very small val- 
ues of J/U, ferromagnetic order becomes stable instead 
of (tt, 0)-antiferromagnetism, and the strongest competi- 
tor in the limit J — >■ is qi = (tt, tt) antiferromagnetism. 
At J = and for small 3/4 < [/ < 1.5, both the FM 
and the (7r,7r)-AF state have pronounced orbital order 
corresponding to cf) = tt/ A: One of the orbitals given by 
the linear combinations Qxz) ± \yz))/V^ is almost filled, 
the other contains sa 1.5 electrons, and xy the remaining 
0.5. In the the FM state, {\xz) + \yz))/V2 is the almost 
filled orbital, while {\xz) — \yz))/^/2has higher occupancy 
in the (tt, 7r)-AF state with only slightly higher energy. 
At larger U > 1.5, both xz and yz are almost filled for 



0.35 
0.3 

0.25 
0.2 

0.15 
0.1 

0.05 




{n,0) 
OD 



AO 



• PM 



'SO^ 



FM 



+ + -H|<xx:jie* X )K X 
k 



FM 



(71,0), FO 



FO 



0.5 



1.5 



u 



FIG. 4: (Color online) Qualitative phase diagram in the J/U 
vs. U plane. The shaded area denotes the stability region of 
the realistic (tt, 0) magnetic ordering. The lines are guides to 
the eye, the dashed line approximately indicates the metal- 
insulator transition. The data points were obtained by com- 
paring the energies of the various phases within the mean- 
field approximation described in App. [X] The meaning of 
the symbols is the following: +: (7r,0)-AF state without pro- 
nounced orbital order (orbital disorder, OD); x: (tt, 0)-AF 
state with (tt, tt) alternating orbital (AO) order with = 0; 
empty squares: (tt, 0)-AF state with (0, tt) orbital stripes (SO) 
{4> = 0); *: (7r,0)-AF state with ferro-orbital (FO) order 
(0 — 0); filled squares: FM with FO ordering tendencies, 
cj) = n/4. This FO order is weaker for small U and larger 
J « 0.2. The filled circles denote parameters that do not sup- 
port magnetically ordered states. For small J, some FO order 
with (/) = 7r/4 is found, similar to the FM phase. The empty 
circles at large U and small J denote similar states without 
magnetic ordering, but with extreme orbital order, where the 
xy orbital is (almost) empty, while xz and yz are (almost) 
filled. 



J — )■ 0, xj/ is nearly empty, and there are hardly any 
unpaired spins that could support magnetically ordered 
phases. Notice that such small ratios of J/U are not ex- 
pected to be realistic for the pnictides, where the onsite 
Coulomb repulsion U is strongly screened)^ With grow- 
ing J/U, the magnetic ordering vector switches to (7r,0) 
ai J/U « 1/5 for small U, and even smaller J/U for 
larger U. 

For U < 2, the (7r,0) magnetic order remains stable 
for aU values of J/U > 1/5 up to J/U = 1/3. At this 
ratio, the Hund's coupling is such that the inter-orbital 
repulsion felt by two electrons in different orbitals, but on 
the same site, vanishes, and we therefore did not consider 
J/U > 1/3. 



C. Evolution of the ground state as a function of U 

for J = U/4 

In this section, the ground state properties at a fixed 
ratio J/U = 1/4 will be discussed. Around this value 
of J/U is where we have found the realistic AF order 
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FIG. 5: (Color online) (a) Orbital magnetization and (b) 
occupation number as a function of the Coulomb repulsion 
strength U, obtained with a mean-field approximation. The 
colors indicate the different phases (for increasing U): uncor- 
related metal, itinerant (tt, 0) antiferromagnet without orbital 
order, itinerant (tt, 0) antiferromagnet with alternating orbital 



order [small white window, spin-orbital order as in Fig. 2(h) 



and a ferro-orbitally-ordered (tt, 0) antiferromagnetic insula- 
tor [spin-orbital order as in Fig. |2(a)| . Hopping parameters 
are from Tab.]!] and J — C//4. For the phase with alternating 
orbital order, the thin lines show (a) m ± g and (b) 2n ± p. 



with ordering momentum qi = (tt, 0) for all values of 
U > Uci- Figure 5(a) shows how the staggered magneti- 



zation with ordering momentum (tt, 0) increases with the 
Coulomb repulsion U . As previously found for two- and 
four-orbital models^ intermediate U leads to an antifer- 
romagnetic metal. The system remains nonmagnetic for 



small U up to Uci ~ 0.6. 



For U > Uci, the spin ( 7r,0)- 



ordered magnetic moment starts to grow, see Fig. 
but the band structure remains metallic, as can be de- 
duced from the spectral functions presented in Fig. El cal- 
culated for several values of the Coulomb repulsion in this 
regime. Note that the spin-(7r,0) AF order triggered by 
U introduces gaps and magnetically-induced "shadow" 
bands)^ for comparison, the uncorrelated A(k, oj) is in- 
cluded in Fig. ini Since several bands are involved, the 
gaps are not necessarily located at the chemical poten- 
tial, as it has been discussed for the two- and four-orbital 
models in Ref. [25l . and the system remains metallic. As 
it can be seen in Fig. [51 the overall features of the spec- 
tral density in this regime remain similar to those of the 
noninteracting limit, with the bandwidth being slightly 
reduced with increasing U. However, the onset of AF 
order docs affect some details of j4(k, w), especially low- 




1,0) (0,7t) (7t,E) (0,0) 





FIG. 6: (Color online) Spectral density yl(k,aj) for the anti- 
ferromagnet orbital-disordered metallic phase at (a) U = 0.7, 
(b) U = 0.9, and (c) U = 1.1. The BZ is for the one-Fe unit 
cell, and J = 11/4 was used. The uncorrelated A{\i,u}) for 
U — is included as solid lines for comparison. 



energy features at the chemical potential, where one of 
the hole pockets disappears and additional pockets arise. 

The Fermi surface for U = 0.7, where the Coulomb re- 
pulsion is just barely strong enough to induce (tt, 0) an- 
tiferromagnetism, is shown in Figs. [7l[a,b). More specif- 
ically. Fig. 7(a) shows the Fermi surface in the extended 
BZ for spin stripes running along the y-direction, i.e., 
for the ordering vector (tt, 0). While the electron pocket 



at (0, tt) is hardly affected, the pocket at (tt, 0) has al- 
most disappeared. Of the two hole pockets, the inner 
one has also disappeared for momenta (0, ky), because a 
gap has here developed at the chemical potential ^, see 
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Fig. 6(a) For momenta (^3^,0), in contrast, the gap in 
the outer pocket Hes below ^, and the band consequently 
forms a very small electron pocket. This result is in qual- 
itative agreement with the unconventional electronic re- 
construction observed with ARPES in (Ba, Sr)Fe2As2i^ 
Figure p7(b) [ shows the superposition of the Fermi surfaces 
obtained for the two equivalent ordering vectors (tt, 0) 
and (0, tt) in the reduced BZ corresponding to the two- 
Fe unit cell. If U is increased to U = 0.9, the gap in the 
outer hole pocket along (fca;,0) increases and pushes the 
outer band above the chemical potential; the small elec- 
tron pockets seen for U = 0.7 in Fig. 7(a) consequently 
disappear, and only one hole pocket remains around F, 
see Fig. 7(c) The (0,7r) electron pocket remains un- 



affected, but at (7r,0), a hole-like shadow pocket with 
very low spectral weight has replaced the original elec- 
tron pocket. The band that formed the vanished electron 
pocket at f7 = has been deformed strongly enough to 
create a small /loZe-like pocket at « (7r/2, 0). As it can be 
seen in Fig. |7(d)[ this hole pocket touches the (tt, 0) elec- 
tron pocket once the results for ordering vectors (tt, 0) 
and (0, tt) are combined. As U continues to increase 
within the magnetic metallic phase no further qualita- 



tive changes are observed, as it can be seen in Fig. 6(c) 
and Fig. 7(e) and |7(f)| where the spectral functions and 
the FS are shown for U = 1.1. 

The average electronic occupation numbers for the 
three orbitals, shown in Fig. 5(b)[ are not significantly 
affected by the onset of antiferromagnctism. We be- 
lieve that the small difference in electronic population 
observed is driven by the different orbital magnetization 
(see Fig. 5(a) ) and is due to the orbital anisotropy relative 
to the direction of the magnetic (tt, 0) stripes. Note that 
the difference between m 
than the difference between n. 



xz and niyz in Fig. 5(a) is larger 



and Uyz in Fig. |5(b) 



dicating that q is more important than p in Eq. (pfT In 
addition, the behavior of the spectral functions appears 
to be dominated by the magnetization in this phase. 

When a second critical coupling Uc^ ss 1.23 is reached, 
the system develops orbital order with an ordering mo- 
mentum (7r,7r), different from the magnetic ordering 
vector (tt, 0); the spin-orbital order realized is the one 
schematically depicted in Fig. 2(h) (the order occurs be- 
tween the xz and yz orbitals, i.e., = 0.) The system re- 
mains a metal through this second transition as well, but 



the spectral density is profoundly affected, see Fig. 8(a) 
The original hole and electron pockets around F and M 
completely disappear and only correlation-induced pock- 
ets remain: the hole pocket around k = (7r/2,0) is mir- 
rored at k « (7r/2, tt), and four more small pockets can be 
seen in the FS (not shown) away from the high-symmetry 
directions plotted in Fig. 8(a) Different states with dif- 



ferent orbital ordering patterns have only slightly higher 
energies in this regime. In contrast, phases with different 
magnetic ordering have significantly higher energy, sug- 
gesting that the (tt, 0) AF "stripes" may be more robust 
than the alternating orbital order. 

If U is further increased, a metal-insulator transition 
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FIG. 7: (Color online) Fermi surface in the orbital-disordered 
spin-antiferromagnetic metallic phase with (a,b) U = 0.7, 
(c,d) U = 0.9, and (e,f) U = 1.1. (a,c,e) show the unfolded 
BZ containing one Fe, for the antiferromagnetic ordering vec- 
tor q = (tt, 0). (b,d,f) depict the superposition of the FSs for 
q = {n, 0) and q = (0, tt) in the (rotated) folded BZ corre- 
sponding to two Fe atoms. The ratio J = U/4: was used. The 
color scale is the same as in Fig. |6l 



finally occurs at a third critical Uc^ ~ 1.43. At this point, 
the orbital order changes: as can be concluded from the 
orbital densities shown in Fig. 5(b) the system develops 
ferro-orbital order. The spin-(7r, 0) antiferromagnctism 
persists, and the spectral density in Fig. |8(b)| has a full 
gap. The ferro-orbital spin-(7r, 0) order in this insula- 

With 



tor is the one depicted schematically in Fig. |2(a) 



growing U, the staggered magnetization converges to its 



maximal possible value 2/iBohr, as shown in Fig. 5(a) 

The direction of the magnetic stripes determines which 
of the two degenerate xz and yz orbitals is (almost) dou- 
bly occupied in the FO order realized at large U. For 
ordering vector (7r,0) ((0, tt)) it is the yz (xz) orbital. 
This can be understood by considering the interorbital 
hopping ti between the xy orbital and the yz{xz) orbital 
along the y-(a;-)direction, which has to be large in or- 
der to ensure xz/yz character for the hole pockets, see 
Sec. Ill Bl For spin stripes along the j/-direction, i.e. or- 
dering vector (7r,0), bonds in the x-direction are anti- 
ferromagnetic and the electrons can better take advan- 
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tage of the AF superexchange if the two orbitals con- 
nected by ^7 are both singly occupied (see the schematic 
iUustration in Fig. The remaining yz orbital then 
has to be doubly occupied. This does not cost any AF 
superexchange energy, because its connection to the xy 
orbital via ty lies along the spin-aligned j/-dircction where 
such AF superexchange would in any case not occur. In 
fact, the additional electron in the yz orbital has the op- 
posite spin from the majority spin of the stripes and can, 
thus, gain some kinetic energy by hopping via t-^ to the 
xy orbital along the FM y-direction. Thus, in this regime 
of large U the ground state corresponds to the cartoon 
shown in Fig. 2(a) if the magnetic order is (7r,0). 

Summarizing, our mean-field calculations indicate the 
existence of four distinct phases that are stabilized with 
growing Coulomb repulsion U: (i) a disordered, param- 
agnetic phase for U < Uci, (ii) a metallic phase with 
(tt, 0) or (0, tt) magnetic order for [7^ < U < Uc^ , (hi) 
a metallic magnetic phase for < U < Uc^ with alter- 
nating orbital order with ordering vector (tt, tt), and (iv) 
a fcrro-orbitally ordered insulator with spin-(7r,0) mag- 
netic order for U > Uc^, where the yz [xz] orbital has 
larger electronic occupation for magnetic ordering vector 
(^,0) [(0,7r)]. 



(7t,0) (0,Jl) (7:,7l) 



FIG. 8: (Color online) Spectral density yl(k,aj) for (a) the 
orbitally-ordered spin-(7r,0) antiferromagnetic metallic phase 
at [7 = 1.36, and (b) the orbitally polarized spin-(7r,0) anti- 
ferromagnetic insulator at (7 = 1.6. The ratio J — 11/4 was 
used, and the unfolded BZ is for the one-Fe unit cell. 



t7[ 




FIG. 9: (Color online) Magnetic order and orbital occupa- 
tion at large U shown for four sites. The magnetic ordering 
vector is (7r,0), i.e, the spin stripes run along the y direc- 
tion. For each site, the xy, xz, and yz orbitals are shown: 
the xy is the one below the other two, and the yz is dou- 
bly occupied. Dashed (continuous) lines indicate inter-orbital 
hopping t7 connecting the xy orbital to xz (yz) along the x- 
(y-) direction. 



IV. PAIRING OPERATORS IN A 
THREE-ORBITAL MODEL FOR PNICTIDES 

In this section, the spin-singlet pairing operators that 
are allowed by the lattice and orbital symmetries in the 
three-orbital model for LaOFeAs will be constructed. 
This classification of operators has previously been made 
for the two-orbital modeli^"— An approach similar to 
Ref. [l^ will be followed. To achieve this goal, the three- 
orbital tight-binding portion of the Hamiltonian, Htb, 
presented in Eq. ^ will be rewritten in terms of the 
3x3 matrices Xi which correspond to the eight Gell'mann 
matricesSS. for the cases i = 1 to 8, while Aq is the 3x3 
identity (see App. |B]). Then, H^b becomes 



ifTB(k) = ^a>^a*k„ 



(14) 



where ^ = (4,(k), eft Jk), 4^(k))^ and 

= fkAo + (^kAs + 7kAi + a'^^ X5 + al^^ Xj + hi^Xs, (15) 



with 



ck = (T" + r22 + r33)/3 

2 

= ^(ii + <2 + t5)(cosfc^ + cosfc.y) 

4 Axy 

+ -(2*3 + te) cos cos ky - fi+ 



(16) 



(17) 



Sk = (r" - r22)/2 = -{h - f,)(cos/.v - COS fey), 

(18) 

= T^^ = 4t4sinfca;Sinfcy, (19) 
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TABLE III: Symmetry properties of the several terms in the 
Htb Hamiltonian of the three-orbital model. 



Term 


IR 




Ala 




Big 


7k 


B2g 




Ea 




Alg 



TABLE IV: Symmetry properties of Gell'mann matrices for 
the orbital assignment that defines the proposed three-orbital 
model. 



Matrix IR 



Ao 


Alg 


Ai 


B2g 


Aa 


A2g 


A3 


Big 



(A4,A6) Eg 

(As, At) Eg 

As Alg 



TABLE V: Form factors /(k) for pairs up to distance (1,1) 
classified according to their symmetry under operations. 



# 


./(k) 


IR 


1 


1 




2 


cos kx + cos ky 


Alg 


3 


cos kx cos ky 


Alg 


4 


cos kx — cos ky 


Big 


5 


sin kx sin ky 


B2g 


6 


(sin kx, sin ky) 


Ea 


7 


(sin kx cos ky , sin ky cos kx ) 


Eg 



TABLE VI: Product table for the irreducible representations 
of the group D4,h relevant to this work. 





Alg 


A2a 


Big 


B2g 


Eg 


Alg 


Alg 


A2g 


Big 


B2g 


Eg 


A2g 


A2g 


Alg 


B2g 


Big 


Eg 


Big 


Big 


B2g 


Alg 


A2g 


Eg 


B2g 


B2g 


Big 


A2g 


Alg 


Eg 


Eg 


Eg 


Eg 


Eg 


Eg 


Alg + A2g + Big + B2g 



^(1) ^ yis/j = _2t7sinfc^ - 4t8sinfc:i, cosfcj^, (20) 



J2) _ ^23 / ■ _ 



T ji = —2t-j sin ky — 4tg sin ky cos k^, (21) 



and 



T + T 



^33 

7f 



(22) 



= — p(tl + <2 — 2t5)(cOS kx + COS ky) 

v3 

-I- —^{ts - te)coskx cosky - —jM-. 

It can be shown that each element in Eqs. (|15ll22p 
transforms according to one irreducible representation of 
the group corresponding to the Fe lattice. The clas- 
sification is given in Tab. IIIIl 

Since the Hamiltonian has to transform according to 
Alg, the GcU'mann matrices in the orbital basis here cho- 
sen transform as indicated in Tab. IIVI 

In multiorbital systems the general form of a spin- 
singlet pairing operator is given hy^ 



At(k) 



/(k)(A0a,/3(4 



7t 



(23) 



where a sum over repeated indices is implied; the oper- 
ators <ij(. Q g. have been defined in the previous sections 
and /(k) is the form factor that transforms according 
to one of the irreducible representations of the crystal's 
symmetry group. Although /(k) may, in general, have 
a very complicated form, a short pair-coherence length 
requires the two electrons that form the pair to be very 
close to each other. Consequently, for simplicity we focus 
on nearest and diagonal next-nearest neighbors, and form 
factors that arc allowed in a lattice with D^h symmetry. 
The momentum dependent expression, as well as the ir- 
reducible representation according to which each form 



factor transforms, are given in Tab. |Vl Note that if the 
pairing mechanism is non-BCS and if the Coulomb repul- 
sion is strong on-site pairing, then f(k)=l corresponding 
to onsite pairing is an unlikely factor. 



A. Intraorbital Pairing 

The previous discussion shows that the symmetry of 
the pairing operator will be exclusively determined by the 
symmetry of /(k) only if Ai transforms according to Aig. 
Table IIVI indicates that this is the case for pairing oper- 
ators constructed by using Ao or Ag in Eq. ([23)) . These 
two matrices are diagonal, which means that such pair- 
ing operators define intraorbital pairings. For intraor- 
bital pairing, with a symmetry fully determined by the 
spatial form factor, the basis functions arc then given by: 
#I:/(k)Ao or #11: /(k)A8. 

For #1, the superconducting order parameter (OP) will 
be the same for the three orbitals while #11 allows the OP 
for xy to be different than for xz and yz which need, by 
symmetry, to have OPs that can only differ by a relative 
sign. Thus, the addition of a third orbital may allow the 
possibility of different superconducting gaps in the band 
representation, reminiscent of the two gaps in MgB2<^i^ 

When any of the remaining seven matrices Ai appear in 
Eq. (|23p . the symmetry of the pairing operator is given 
by the irreducible representation of D^h resulting from 
the product of the symmetry of the form factor and the 
symmetry of the orbital component?^ according to the 
product table given in Tab. Ell For Ai = A3, the basis 
function is given by #111: /(k)A3. The pairing is still 
intraorbital but since A3 transforms according to Big the 
symmetry of the operator will be Big if /(k) transforms 
according to Aig, etc. Note that this pairing operator 
does not involve the xy orbital and it has already been 
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presented in the context of the two-orbital modeli^ How- 
ever, since the orbital composition of the bands is not the 
same as for the two-orbital model, it will be important 
to determine whether the gap structure of this pairing 
operator has changed. 

B. Interorbital Pairing 

The remaining six A; matrices lead to interorbital pair- 
ing. Note that Ai and A2 do not involve the orbital 
xy and the pairing operators that they generate have 
already been discussed in the two-orbital model.— We 
are interested in the spin-singlet pairing operator for or- 
bitals xz/yz that has a basis 7^IV:/(k)Ai. This opera- 
tor, with /(k) = coskx + cosky, has been found to be 
favored for intermediate values of the Coulomb repul- 
sion U in numerical calculations of the two-orbital model 
for pnictides<^i2i The addition of the xy orbital leads 
to the possibility of new interorbital pairing operators, 
i.e., pairing between electrons in the orbitals xz and yz 
with electrons in the xy orbital. Thus, now the focus will 
be on the interorbital spin-singlet pairing operators that 
result from the addition of xy. 

The interorbital case becomes very interesting because 
we need to combine {xz,yz) that transform as the two- 
dimensional representation Eg with xy that transforms as 
B2g; thus, the product transforms as Eg. The A^ matrices 
that can appear in this intraorbital pairing are (A4,A6) 
or (A5,A7). Since the focus here is on pairing operators 
that are spin singlets, it will be required that the oper- 
ator is even under orbital exchange. Thus, only (A4,A6) 
will be considered since the other two matrices (A5,A7) 
that transform according to Eg will produce operators 
odd under orbital exchange. Let us further restrict the 
analysis to the case of pairing operators that transform 
according to one dimensional representations of the point 
group because we assume that the ground state is non de- 
generate. Thus, the only spatial form factors /(k) that 
we should consider must transform according to Eg . This 
leaves us with 

/(k) = (sinfc:r,sinfcy), (24) 

for nearest neighbor pairs and 

/(k) = (sin fc^ cos fcj^, cos fca; sin fcj,), (25) 

for diagonal pairs. Since the direct product of two Eg 
representations is EgX Eg ~ Aig + A2g + Big + B2g, pair- 
ing operators transforming according to four irreducible 
representations will be obtained. The basis for the new 
pairing operators, labeled Vi, are presented in Tab. IVIll 

C. Band Representation 

To obtain the gap structure of the pairing operators #1 
to #IV (shown in Table |VlT| the Bogoliubov-de Gennes 



TABLE VIL Properties of pairing operators in the three- 
orbital model. / indicates the symmetry of /(k). 



No. 


IR 


Basis 


Gap 


I 


f 




/(k)Ao 


Full or Nodal 


II 


f 




/(k)A8 


Full or Nodal 


III 






.f(k)A3 


Nodal 


IV 






/(k)Ai 


Nodal 


Ya 


Ala 




sin fca; A4 + sin fcj, Ag 


Nodal 


Vb 


Ala 


A4 


sin cos ky -f Ag cos k^ sin ky 


Nodal 


Vc 


Big 




sinfca;A4 — sinfcyAe 


Nodal 


Vd 


Big 


A4 


sin kx cos ky — Ae cos k^ sin ky 


Nodal 


Ve 


Alg 




sin fca; Ae + sin fcj, A4 


Nodal 


V/ 


A2g 


A4 


cos kx sin ky -f Ag sin kx cos ky 


Nodal 


V, 


Big 




sinfc^Ag — sinfcj,A4 


Nodal 


Vh 


Big 


A4 


cos kx sin ky — Ag sin kx cos ky 


Nodal 



Hamiltonian (BdG) is constructed and it is given by 

i/BdG = E*W^*k, (26) 

k 

with the definitions 

*k = ("^k.xz.t' ^k,i/2,t' ^k,2;y,t' 

d-]i^xz,i, d-k^yz,l, d-k^xy,i), (27) 

and 

„MF_ /^i^TB(k) P(k) \ 

[ P^k) -ifTB(k) J ' ^^^> 

where each element represents a 3 x 3 block with Htb (k) 
given by Eq. (|4]) and 

P(k)„,^ = y/(k)(A,)a,/j, (29) 

with i = 0, 8, 3, and 1 for pairing #1, #11, #111, and 
#IV, respectively. V is the magnitude of the OP given by 
the product of the pairing attraction Vq and a mean-field 
parameter A that should be obtained from minimization 
of the total energyi^ For pairing #Vi the basis listed in 
Table IVlIl should be used instead of /(k)Ai. 

Up to this point we have worked using the orbital 
representation because this basis renders it straightfor- 
ward to obtain the form of the Hamiltonian, as well as 
the pairing operators allowed by the symmetry of the 
lattice and orbitals. However, the experimentally ob- 
served superconducting gaps occur at the FS determined 
by the bands that result from the hybridization of the 
orbitals. For this reason, it is convenient to express 
Eq. pS]) in the band representation. HtbO^) can be ex- 
pressed in the band representation via the transformation 
i?Band(k) = J7t(k)7JTB(k)f7(k), where U(k) is the uni- 
tary change of basis matrix and J7^(k) is the transpose 
conjugate of J7(k). Since U is unitary it is known that 
for each value of k, J2iiUh3)*Ui,k = Y.i{Uj,i)*Uk,i = 5j,fe. 
Then, -ff^F ~ HmfG where G is the 6x6 unitary ma- 
trix composed of two 3x3 blocks given by U. Then, 

rr/MF _ 

/i?Band(k) Pij(k) \ 
-[ P^(k) -i/Ba„d(k)j' 
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with 

PB(k) = U-\k)P{k)U{k). (31) 

A standard assumption in superconducting multiband 
systems is that the pairing interaction should be purely 
intraband, meaning that PbO^) is diagonal. Thus, let 
us explore what kind of purely intraband pairing opera- 
tors are allowed by the symmetry properties of the three- 
orbital model for LaOFeAs. In the band representation, 
the most general BdG matrix with purely intraband pair- 
ing is given by 



Ai(k) 








Ai(k) 










e2(k) 








A2(k) 











e3(k) 








A3(k) 


Alik) 








-ei(k) 



















-e2(k) 





V 





A5(k) 








-e3(k)/ 



(32) 

where ei(k) are the eigenvalues of iJTB(k) and Ai(k) 
denotes the band and momentum dependent pairing 
interactions!^ As it can be deduced from the properties 
of the unitary change of basis matrix U, if all three bands 
have the same pairing interaction, i.e. Ai(k) — A2(k) — 
A3(k) = A(k), then the matrix P{k) in the orbital rep- 
resentation will also be diagonal. In this case, the pairing 
operator is given by Eq. with an arbitrary /(k) and 
\ = ^0, i-e., the pairing operator is intraorbital and the 
OP is the same for the three orbitals. This corresponds 
to pairing operator ^I which describes a pairing interac- 
tion that is the same for each of the three orbitals. How- 
ever, symmetry only requires that the orbitals xz and yz 
must have the same OP, while xy can have a different 
one. Thus, there does not seem to be a reason to assume 
that electrons in the many bands that determine the FS 
should be affected by the same pairing interactions. In 
fact, in MgB2 the electron-phonon interaction that pro- 
vides the pairing is stronger on the cr-bands than on the 
TT-bands giving, as a result, two different superconducting 
gaps. Thus, it can be asked whether the symmetry of the 
three-orbital model allows for the possibility of two dif- 
ferent OPs with a pure intraband pairing interaction. If 
it is assumed that in Eq.(132]) Ai = A2 = A = /(k)C and 
A3 = A' = /(k)C", then in the orbital representation 

PikUf3 = f/o,3(k)C/^,3(k)/(k)(C' - C), (33) 

for the off-diagonal elements and 

P(k)„.„ = A + |[/„,3(k)P/(k)(C' - C), (34) 

for the diagonal ones. 

Now let us concentrate on the diagonal part. This has 
to arise from a linear combination of intraorbital pairing 
operators with compatible symmetries. There are two 
possibilities: 

P(k)a,a = /(k)[A(Ao)a,a + B{Xs)a.,c.], (35) 



or 

P(k)„,„ =I?/(k)(A3)„,a, (36) 

where A, B, and D are independent of momentum. It 
can be shown that Eq. (|35|) requires |C/i3p = |C^23p while 
Eq. dSni) requires |t/i3p — — IC/23P, which are not satisfied 
by the elements of the matrix U determining the change 
of basis. This means that any purely intraband pairing 
interaction allowed by the symmetry of the three-orbital 
model should be the same for the three bands. On the 
other hand, if we had a case in which |?7i3p = IC/23P = 
0, which means that one of the three orbitals does not 
hybridize with the other two, it would be possible to have 
a system with two different gaps. Note that this is the 
situation for MgB2 in which the z orbital that forms the 
TT band does not hybridize with the x and y orbitals that 
constitute the a band. 

Summarizing, it has been found that independent gaps 
in different Fermi surfaces cannot arise if the hybridiza- 
tion among all the orbitals is strong and the pairing in- 
teraction is purely intraband. 

D. The s± Pairing Operator 

The next issue to be considered is whether the orbital 
and lattice symmetries allow for the possibility of the of- 
ten discussed s± pairing scenario. ARPES experiments 
indicate the existence of two hole-pockets around F. The 
interior pocket, which is almost nested with the electron 
pockets with a nesting vector q = (tt, 0) or (0,7r), devel- 
ops a constant gap Ah which has the same magnitude 
than the gap on the electron pockets Ag. In addition, 
they find a smaller gap A'^ « A;i/2 on the exterior hole 
pockct41 The ARPES resuhs can be interpreted in two 
different ways in the context of a three-orbital model: 

(i) Assume that the inner hole pocket observed in 
ARPES corresponds to two almost degenerate FS that 
cannot be resolved, and assign the external hole pocket 
to a band that arises when extra orbitals are added. This 
is the same assumption made in the two-orbital model for 
which it was shown in Ref. [l^ that the s± pairing state 
is compatible with the lattice and orbital symmetries. 

Under this assumption, in the three-orbital model the 
s± pairing state corresponds to our pairing operator #1 
with /(k) = cosfc^: cos ky, which in the band representa- 
tion leads to a purely intraband pairing attraction given 
by Ai{k) ~ V cos kx cos ky for each of the three-bands. 
For hole pockets almost degenerate with each other, the 
gap in both bands will be the same and there will be a 
sign difference with the gap at the electron pockets whose 
Fermi momentum differs from those of the hole pockets 
by (0,^) or (^,0). 

(ii) Assume that the inner and outer hole pockets ob- 
served by ARPES are described by the two hole pock- 
ets in the three-orbital model. This would force us to 
request that, for example, Ai(k) = — A3(k-|-q) where 
q = (7r,0) or (0,7r) and A2(k) is independent. Then, let 
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us assume that Ai(k) = A3(k) = cos fc^^ cos fc^Ao and 
A2(k) = A(k). Let us concentrate on the T-X direc- 
tion. Along this direction, 7k and a^"^ vanish, meaning 
that there is no hybridization among the three orbitals. 
From Fig. [1] we observe that each of the two hole FS re- 
sults from the crossing of xz and yz, while the electron 
FS has pure xy character. From the orbital symmetry, 
then, it is deduced that the only reason for having dif- 
ferent gaps at the two hole-like FS would be a strong 
momentum dependence of the gap since symmetry en- 
forces |A2;z(k)| = |Ay2(k)|; in addition, the gap at the 
electron pockets does not need to be related to the gap in 
the hole pockets, unless the pairing operator contains Aq. 
Thus, we observe that the s± pairing operator could be 
supported under this assumption if it is given by Eq. ([23]) 
with /(k) = cos kx cos ky with = Aq and the additional 
condition that if k'^'^'' represents the Fermi momentum 
of the internal hole pocket and k^^f^' = k^'^'' -I- (5 is the 
Fermi momentum of the external hole pocket it is nec- 
essary that /(k^'^'')//(k^'^''' ) w 2 which would require 
fine tuning of the parameters. It would also be expected 
that in this scenario the ratio A/j/ /Ah should not be 1/2 
for all materials. 

Thus, it is concluded that the s± pairing could be sup- 
ported by a three-orbital model. It corresponds to pairs 
of electrons in the same orbital at distance one along 
the diagonals of the square lattice, i.e., on next- nearest 
neighbor sites, with the same pairing potential for all 
three orbitals. Then, if experiments show that s± is in- 
deed the correct pairing operator, it will remain to be 
understood why the pairing interaction does not appear 
to depend on the symmetry of each different orbital or, 
equivalently, why it is the same for electrons in different 
bands. This should be contrasted with the case of MgB2 
in which the strength of the eleetron-phonon coupling 
that leads to pairing is stronger on the cr-band FS than 
in the 7r-band FS. 

Since the pairing mechanism for the pnietides is not 
known and the s± pairing state is just one of many pro- 
posed states, our discussion will continue by analyzing 
the other new pairing states that are allowed by symme- 
try when the xy orbital is considered. 



E. Properties of the Pairing Operators 

In the previous subsection, it was shown that only pair- 
ing #1 can lead to a purely intraband pairing interaction 
in the context of the three-orbital model, and that the 
widely proposed s± pairing state indeed belongs to the 
class represented by pairing ^I. On the other hand, given 
the complexity of the problem, only numerical simula- 
tions can clarify whether the three-orbital model becomes 
superconducting upon doping and what is the symmetry 
of the dominant pairing state. Since numerical results are 
not yet available, here the properties of the other possible 
spin singlet pairing operators will be discussed. 

Let us start with pairing #11, i.e. the intraorbital 



pairing operator containing As, that allows a different 
pairing strength for the xy orbital. This pairing op- 
erator is not purely intraband. This fact can be eas- 
ily deduced from the properties of the unitary change 
of basis matrices. In this ease, P(\i)a.fi = CaSafj with 
Ci = 6*2 = l//(k)/\/3 corresponding to orbitals xz and 
yz, and C3 = — 2F/(k)/\/3 for orbital xy. Then, 

PsikUfi = ^aC/*„(k)C/,,^(k), (37) 

i 

which does not vanish for all values of k for a 7^ /3, thus 
indicating the existence of interband pairing terms. Sim- 
ilar calculations for all the pairing operators presented in 
Table IVlIl show nonvanishing interband pairing terms. 

Then, it is concluded that starting from the orbital 
representation, the only way to obtain pure intraband 
pairing in the band representation is by considering a 
pairing interaction that affects equally all the orbitals in- 
volved, producing equal gaps in all the orbitals and/or 
bands with symmetry determined by the spatial form 
factor. This shows that the requirement of purely intra- 
band pairing induces a strong constraint regarding the 
coupling of the electrons in the different orbitals with 
the source of the pairing attraction. On the other hand, 
if the requirement is relaxed, interband pairing occurs at 
least in some regions of the BZ.— It was verified that this 
is the case for the remaining operators #11, #111, #IV, 
and #Vi. It has also been observed, by monitoring the 
eigenvalues of -ffedc for operators #111 and #IV, that 
there is a nodal structure in the superconducting gap for 
all the values of /(k) shown in Tab. |Vl while operator 
#11 becomes nodeless for /(k) = cosfc2;COsfcj, or 1 at a 
finite value of V. In addition, some linear combinations 
of pairing #1 and #11 with /(k) = cosk^ cos ky or 1 are 
also nodeless for all finite values of V but they lead to in- 
terband pairing interactions. These nodeless states that 
we call Sib will be discussed in Sec lIV E"2] 



1. Pairing with Pseudocrystal Momentum Q = (tt, tt) 

In Sec. IIIBl it was explained that although the three- 
orbital Hamiltonian retains the two-iron unit cell of the 
original FeAs planes, it is possible to express it in terms 
of three orbitals in the space of pseudocrystal momen- 
tum k defined in the extended Brillouin zone corre- 
sponding to a base with one single Fe atom per unit 
cell. In terms of the real momentum, the Hamilto- 
nian consists of two 3x3 blocks Hi (k) and H2 (k) with 
i/i(k) ^ i?TB(k) = i72(k + Q); thus, these two blocks 
provide the same eigenvalues in the unfolded Brillouin 
zone, but both blocks need to be considered if the actual 
reduced BZ is used. This means that in the reduced BZ 
the bands arise as combinations of six orbitals labeled by 
the orbital index a ~ 1. 2, or 3, and the Hamiltonian 
block index i = 1 or 2. Then, to consider all the possible 
interorbital pairing operators it is important to include 
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pairs formed by electrons in orbitals in the two different 
blocks. In the extended BZ, this is equivalent to con- 
sidering pairs with both pseudocrystal momentum and 
Q. Note that Cooper pairs with pseudocrystal momen- 
tum Q still have zero center-of-mass momentum. Exact 
diagonalization studies of the two-orbital model^i^ did 
not favor such operators, and it is possible that this kind 
of pairing docs not occur in the three-orbital model ei- 
ther. However, since symmetry allows such a possibility, 
pairing operators with nonzero pseudocrystal momentum 
will here be discussed for completeness. 

The generalized Bogoliubov-de Gennes matrix 
that allows us to consider interorbital pairs with pseu- 
docrystal momentum Q is given by: 



MF 



/i?TB(k) P(k) \ 

-iJTB(k) F(k + Q) 
pt(k + Q)ifTB(k + Q) 
\P^(k) -HTB{k + Q)J 

(38) 



where P(k) has the form given in Eq. ([29l) . 

By finding the eigenvalues of , the structure of the 
gap of the possible pairing operators with pseudocrystal 
momentum Q can be obtained. Our analysis shows that 
all the pairing operators with pseudocrystal momentum 
Q lead to inter and intraband pairing in the band rep- 
resentation and nodes on the FS for small V. We have 
observed that a nodeless gap for pairs with pseudocrystal 
momentum Q develops at a finite value of V for operators 
and #11 with /(k) ~ 1 or cosfc^^cosfej^ in a manner 
characteristic of systems with interband pairing.— 



2. Spectral Functions 

It is straightforward to calculate the spectral func- 
tions A(k, oj) for all the pairing operators presented in 
this manuscript. However, due to their large number, 
we will concentrate on (i) pairing operator #1 with 
/(k) = cosfcj;Cosfcy, i.e., the s± pairing operator, (ii) 
a linear combination of pairing operators #1 and #11 
with /(k) = cosfc^ cosfcy, that we will call the sjb pair- 
ing operator; (Hi) pairing operator #IV with /(k) = 
cos kx + cos ky, i.e., the i?2g pairing operator, favored by 
numerical calculations in the magnetic metallic regime 
of the two-orbital model, which will be called i?2g, and 
(iv) a, linear combination of pairing operator #IV with 
/(k) = cosfc^, -I- cosky and pairing operator with 
pseudocrystal momentum Q, which is the natural exten- 
sion to three orbitals of B2g and will be called i?2g* ■ 

In Fig. rTOf a) the spectral functions A(k, oj) along high 
symmetry directions in the reduced Brillouin zone are 
shown for the three-orbital Hamiltonian with V = 0, i.e., 
without pairing, in order to illustrate the changes induced 
by the various pairing interactions considered here. Note 
that these results correspond to the system with two Fe 
atoms per unit cell, which leads to the six bands seen 
in Fig. [TUfa). The results for the s± pairing operator 
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FIG. 10: (Color online) The intensity of the points represents 
the values of the spectral function A(k, lj) for the three-orbital 
model with pairing interaction (a) V = 0; (b) V — 0.2, for 
the s± pairing operator given in the text. 



with intensity V = 0.2 are presented in Fig. [TOT bl. It 
can be observed that a gap opens at the FS and shadow 
Bogoliubov bands, which should be visible in ARPES 
experiments, appear. Numerically, we have verified that 
no nodes occur anywhere in the BZ. Note that the gap 
is momentum dependent because /(k) ~ cos k^ cosky. 
This means that the ratio between the gaps at the FS 
is determined by |/(ki)//(kj)|, where i and j can take 
the values 1, 2, 3 corresponding to the three bands that 
determine the FS, i.e., 1 (2) for the interior (exterior) 
hole pocket, and 3 for the electron pockets. This creates 
a constraint on how different these gaps can be if s± 
represents the actual pairing symmetry of the pnictides. 

As previously mentioned, we can define a pairing op- 
erator Sib by combining operators #1 and #11 with 
/(k) = cos kx cosky so that the orbital part of the ba- 
sis is given by ^Aq + BXg, where A and B are constants. 
This pairing operator is diagonal in the orbital represen- 
tation, and transforms according to Aig thus, it has 5" 
symmetry; it also has intra and interband terms in the 
band representation. This is why this operator is called 
sjB- For a robust range of values of A and B, a nodeless 
gap opens on all FSs for any finite value of V . For ex- 
ample, we can choose the parameters in such a way that 
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FIG. 11: (Color online) (a) The intensity of the points repre- 
sents the values of the spectral function A(k, ui) for the three- 
orbital model with pairing interaction V — 0.2 for the pairing 
state Sib along the indicated high symmetry directions in the 
folded BZ. (b) Ratio R between the gaps for pairing sib and 
pairing s± for V = 0.05 in the unfolded BZ. 



for any given k, the pairings for the three orbitals have 
the same sign. The speetral funetions for y = 0.2 along 
high symmetry directions in momentum space are shown 
in Fig. 11(a) for the special case in which A — 3/2 and 
B = -V3/2. The major difference with the results for 
the s± state [see Fig. rTOT b)] is that the interband pair- 
ing present in s/b opens gaps between the bands away 
from the FS. This is a feature that should be observed in 
ARPES experiments. Also the band spectral functions 
in both cases are very different close to (0, tt) and (tt, 0). 

We have also investigated how the gaps on t he diff er- 
ent FSs differ for pairing s± and s/b. In Fig. 



11(b) 



show the ratio R between the two gaps in the unfolded 
BZ for V = 0.05. It can be seen that on the hole pockets 
R = 1, but an appreciable difference is observed on the 
electron pockets where i? = 2 at the point where the elec- 
tron pocket is entirely formed by xy, and diminishes as 
the hybridization of xy with xz or yz becomes stronger. 
The maximum value of i? is a function of the values of A 
and B in the linear combination that defines sjb- Thus, 



while s± is characterized by gaps with a weak momen- 
tum dependence and with similar magnitudes on the hole 
and electron pockets, the sjb state is characterized by a 
different gap on the electron pockets with stronger mo- 
mentum dependence due to the hybridization. 

Now we focus on the spectral function for operator i?2g 
presented in Fig.fT^a). This is the pairing operator that 
was favored by numerical calculations in the intermedi- 
ate U regime of the two-orbital model and it only pairs 
electrons in orbital xz with electrons in orbital yzM- Al- 
though neither this pairing operator nor the Fermi sur- 
faces defining the two hole pockets involve the xy orbital, 
the results around the hole pockets differ in the two- and 
three-orbital models. In the latter, much lower values 
of the pairing attraction V are sufficient to remove the 
extra nodes found close to the hole-pocket FS along the 
r-X{Y) directions in the two-orbital model. This hap- 
pens because the bands forming the two hole pockets are 
now degenerate at F and the pockets are consequently at 
almost the same momenta of the extended BZ, while they 
were separated by (tt, tt) in the two-band model. As a re- 
sult, a small interorbital pairing can now overcome the 
separation between the two FSs and induce a full gap at 
the hole pockets.— At the electron pockets a third node, 
in addition to the two of the two-orbital model, is found 
along the r-X{Y) direction, where the pocket has purely 
xy character and is, thus, not affected by the operator 
B2g. Thus, this pairing operator would show full gaps on 
the hole pocket FS and nodal gaps on the electron pocket 
FS. 

Figure I12f b) shows the spectral functions for pairing 
operator i?2g*j which is a linear combination of B2g with 
the pairing operator Vg (with crystal momentum Q), i.e. 
next-nearest-neighbor interorbital pairing among elec- 
trons in all three orbitals is allowed. We find that nodes 
occur only at the electron pockets. As V increases, nodes 
at the electron pockets remain only along the r-X(V) 
directions because one of the electron pockets is formed 
by a non-hybridized orbital xy along this direction, and 
the relevant pairing interaction is zero for the Fermi mo- 
mentum. We also investigated A(k, w) for a similarly ex- 
tended Big pairing Bf^^ = {coskx+cosky)X3 — a{coskx — 
cos fcj,)(Ao — V^Xs)/^ (not shown), which is nodal for ex- 
tremely small y, with nodes in the T-M direction on the 
hole pockets, but where the nodes are already lifted for 
finite but small V > 0.01 for many non-zero values of a. 

Summarizing, we have found that among the nearest 
and next-nearest neighbor pairing operators allowed by 
symmetry only the s± pairing operator is purely intra- 
band and produces nodclcss gaps for all values of the pair- 
ing attraction V. Thus, purely intraband pairing interac- 
tions occur only if electrons in each of the three orbitals 
arc subjected to the same pairing attraction, i.e. when 
the identity matrix Aq characterizes the orbital portion 
of the pairing operator. We also found that some linear 
combinations of pairing operators #1 and =ffll produce 
nodeless gaps for any finite V if /(k) = cosk^ cos ky, but 
interband attraction appears in parts of the BZ. Finally, 
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FIG. 12: (Color online) The intensity of the points represent 
the values of the spectral function ^(k, uj) for the three-orbital 
model with the pairing interaction V = 0.2, for the pairing 
operators (a) -629 and (b) B2^^ discussed in the text. 

the intcrorbital pairing operators B2g and i?2g' favored 
by numerical studies in a two-orbital model, present 
nodeless gaps on the hole pockets but nodes appear on 
the electron pockets. 

V. CONCLUSIONS 

In this work, a simple three-orbital Hamiltonian has 
been constructed involving the 3d orbitals xz, yz, and 
xy. These orbitals have the largest weight at the FS of 
the pnictide LaOFeAs, according to LDA calculations. It 
was shown that it is possible to qualitatively reproduce 
the shape of the LDA-FS by fixing the electron filling 
to 4 electrons per Fe. Moreover, two features that have 
been criticized in the two-orbital model have now been 
corrected: both hole pockets now arise from bands degen- 
erate at the F-point, and there is no pocket around AI in 
the extended BZ. In addition, the xy character of a small 
piece of the electron pockets is now properly reproduced. 

Numerical calculations using a small 2x2 lattice show 
a tendency to the development of magnetic (tt, 0)-(0, tt) 
stripes when Coulombic interactions are added, result 
consistent with experimental observations. A mean-field 
analysis confirms this tendency for physically relevant 



values of J /U . As in the case of the two-orbital model, 
an antiferromagnetic metallic phase occurs only at in- 
termediate values of the Coulomb repulsion. At large 
U , the ground state is magnetic, but it is an insulator 
that is also orbitally ordered. Additionally, a metallic, 
magnetic and orbitally ordered phase is encountered just 
before the metal-insulator transition. In the most inter- 
esting regime with a spin-(7r,0) antiferromagnetic metal 
without pronounced orbital order, the bands are similar 
to the uncorrelated ones, but their bandwidth is reduced 
with increasing U . The Fermi surface is also very simi- 
lar to the uncorrelated one but, depending on [/, we find 
small additional electron-like pockets near the original 
hole pockets around F (small U) or hole-like pockets be- 
tween the electron- and hole-pockets (at slightly larger 
U). 

The possible pairing operators that are allowed by 
the symmetry of the lattice and the orbitals have been 
constructed for pairs made of electrons separated by a 
distance up to one diagonal lattice spacing. If on-site 
pairing is disregarded due to the large Coulomb repul- 
sion, it was found that the only purely intraband pair- 
ing operator that has a full gap on the FS is #1 with 
/(k) = coskx cos ky which corresponds to the s± pairing 
operator with a momentum dependent OP that has op- 
posite signs on the hole and electron FSs. This operator 
arises from a purely intraband pairing attraction equal 
for each of the three bands. Note that the pairing op- 
erator #1 is the only one that leads to purely intraband 
pairing interactions. Since this pairing operator is pro- 
portional to the identity matrix Aq both in the orbital 
and the band representations we found that the ratio 
|Ai/Aj| between the gaps in two different FSs can differ 
only by the ratios |/(ki)//(kj)|; then, any experimental 
indication of a different kind of ratio would indicate some 
degree of interband pairing.— Thus, order parameter ra- 
tios predicted by several authors^2i22r2^ with calculations 
based on purely intraband pairing (they allow interband 
hopping of intraband pairs) are not allowed by the sym- 
metry of the lattice and the orbitals. In this regard, our 
calculations seem to indicate that unrelated gaps in dif- 
ferent FSs can occur only in systems in which at least one 
orbital (or a group of orbitals) is not strongly hybridized 
with the remaining ones. 

We found that all the other pairing operators, except 
for 7^1, lead to interband pairing attraction in the band 
representation. In addition, all the pairing operators with 
interband pairing studied here have nodal band struc- 
tures at small V with the exception of pairing operator 
Sib- In this case, the gap on the electron pockets is ex- 
pected to have a stronger variation at different points in 
the BZ that the gap at the hole pockets. Thus, a strong 
indication that s± is the appropriate pairing symmetry 
would be provided by experiments in the pnictides show- 
ing a nodeless gap in all FSs, relatively independent of 
momentum, and with similar values on all FSs. 

Summarizing, we have shown that the addition of a 
third orbital corrects the shortcomings pointed out in 
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the two-orbital model: the two hole pockets now arise 
from bands degenerate at the F point while the elec- 
tron pockets contain a small piece with xy character. 
However, the dependence of the magnetic phases with 
U for the undoped case appears to be similar for three 
and two orbitals except for a magnetic, orbital ordered, 
metallic phase that appears in the three-orbital case. 
In both models it is found that the only pairing oper- 
ator allowed by symmetry with next or diagonal nearest- 
neighbor interactions which is purely intraband and pro- 
duces a nodclcss gap is the s± state. In addition, the only 
change observed in the interorbital B2g pairing state, fa- 
vored by numerical simulations in the two-orbital model, 
is that, at the mean- field level, the addition of the xy 
orbital renders the gap on the hole pockets nodeless for 
much smaller values of the pairing attraction. 
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Appendix A: Mean-Field Equations 

In this appendix, we discuss the mean-field approach 
used here to study the Hamiltonian given by the ki- 
netic energy Eq. ^ and the onsite Coulomb interaction 
Eq. (fTTj) . Depending on the ordering vectors qi and q2, 
listed in Table [TTl for magnetic and orbital order, the 
real-space unit cell contains one, two or four sites: One 
for the ferro-orbital and ferromagnetic case, four if both 
ordering vectors are different from each other and from 
(0, 0), and two in all other cases. Following Refs. [25lf59l . 
we only keep the expectation values of density operators, 
as given in Eqs. (fT2|) and (fT3|) in the Coulomb interac- 
tion Eq. pT|) . Together with Eq. (jl]), this determines 
the mean-field Hamiltonian, see below. We then solve 
the self-consistency equations for the six parameters n^y, 
ruxy, n, m, p, and q in Eqs. (fT2|) and ([T3|) for various 
combinations of ordering momenta, see Tab. |TT1 which 
corresponds to minimizing the total energy. This is done 
for all considered phases and the one with the lowest en- 
ergy is taken to be the stable solution. Depending on the 
size of the unit cell, one to four momenta are coupled by 
the Coulomb interaction. In the following, wc will pro- 
vide the Hamiltonians for several ordering patterns with 
different unit cells. In all cases, the sums run over the 
whole extended BZ corresponding to the one-iron unit 
cell. The calculations were carried out in momentum 
space for up to 400 x 400 k-points. We did not observe 
any pronounced dependence on the number of momenta, 
except for very small lattice sizes. 



1. Ferromagnetic and Ferro-Orbital Order: 
One-Site Unit Cell 

In this case qi = q2 = (0, 0) and 

ffMF(k) = HTB(k) + UJ2 ^A.^.adk,,,. 

+ (2C/'-,/) ^A,,.ad^.,,a 

-f^E ^^A.,.JK,,<r (Ai) 

where = n^y, = m^y for the xy orbital, and = 
n±p/2, = m±q for xz and yz. The sum over k runs 
through the whole BZ, N is the number of lattice sites, 
and the constant C is given by 

C=-UY,nl + U/4Y,K (A2) 
— (2U' — J) ^ n^Tii, + J/4 ^ m^rrii, . 

2. Antiferromagnetic and Ferro-Orbital: Two-Site 
Unit Cell 

For AF order with qi = (7r,7r), (0,7r) or (7r,0), the 
real-space unit cell doubles, and momenta k and k -|- qi 
are coupled by the interaction. 

iJMF(k) = i/TB(k) + UJ2 ^A.^Jk,,^ 

^^A+^,,^,.di,,f.,a (A3) 

k,/i,f7 

Again, = n^y, ~ rrixy for the xy orbital, and 
= n ± p/2, rrifi = m ± q for xz and yz; and the same 
constant Eq. (|A2[) as above. The case of ferromagnetic 
order and alternating orbitals is treated in an analogous 
manner. 



3. Antiferromagnetic and Alternating Orbital 
Order with the Same Ordering Vector: Two-Site 
Unit Cell 

In some phases, both the orbital and the magnetic or- 
der alternate with the same ordering vector q = qi = 
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q2 = (tTjTt), (0, tt) or (7r,0). In this case, the Hamilto- 
nian is given by 

i/MF(k) =i/TB(k) 

+ [(4C/' - 2J)n + Un.,y] ^ dL,,.rfk,.,,. 



[{U + J)m + Jm.^y\ 54+qi,M.^'^k,M,o 



k,o- 

^—xz,yz 



{Um^y + 2Jm) J2 f 4+qi.a.y,a4,^y,a 
k,cr 



k,(T 

fi—xz^yz 

-iU + J)q E ^<,.X,,,. (A4) 



{U-2U'~J)p J2 f4 



k,cr 

aa 



2 k+q2,/i,CT k,/^,cr 



k,cr 

fx—xz,yz 



/J.— a;2,y2: 

- [(C/ + J)m + Jm^y] E f 4+qi,M,a^k,/.,a 

k,<T 

fi—xz,yz 

- {Um^y + 2Jm) E |4+qi,xy, 

k,(T 



The constant C is still given by equation (|A5 



Appendix B: Ai matrices 



([/ — 2J7 — J)p "^'^k+qi,/i,cr'^k,/i,(T + The Ai matrices used in the text are presented here: 



k,(T 

fi—xz^yz 



Here, a = ±1 distinguishes between the xz and yz orbital 
flavors as g docs for the spin. The constant C reads 

- [W - 4:J)n^yn - {4U' - 2J){n^ - p^/A) 

+ Jm^yTn + J{in^ -q^)/2 . (A5) 
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4. Antiferromagnetic and Alternating Orbitals with 
Different Ordering Momenta: Four-Site Unit Cell 

If both orbital occupation and magnetic order alter- 
nate with different ordering momenta, so that q.; = 
(tt, tt), (0,7r) or (tt, 0) with qi 7^ q2, the real-space unit 
cell contains four sites, and consequently all four mo- 
menta k, k -|- qi, k -I- q2, k -I- qi -I- q2 are coupled, but 
apart from this, the Hamiltonian is very similar to the 
previous case: 

i/MF(k) =i/TB(k) (A6) 

+ m' - 2J)n + Un.,y] E 4..,..'^k,.,,. 

k,o- 

+ [Un+{2U' -J){n,,y+n)] ^ 4,^,A.p,. 
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